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Abstract 

The so-called level crossing analysis has been used to investigate the empirical data set. But 
there is a lack of interpretation for what is reflected by the level crossing results. The fractional 
Gaussian noise as a well-defined stochastic series could be a suitable benchmark to make the level 
crossing findings more sense. In this article, we calculated the average frequency of upcrossing 
for a wide range of fractional Gaussian noises from logarithmic (zero Hurst exponent, = 0), 
to Gaussian, H = 1, {{) < H < 1). By introducing the relative change of the total numbers of 
upcrossings for original data with respect to so-called shuffled one, IZ, an empirical function for 
the Hurst exponent versus TZ has been established. Finally to make the concept more obvious, we 
applied this approach to some financial series. 
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I. INTRODUCTION 



We are living in a world in which processes of stochastic type are ubiquitous. Although 
the random values of a stochastic process at different times may be independent random 
variables, in most commonly considered situations they exhibit complicated statistical cor- 
relations. Thus, over past decades several different methods have been introduced to inves- 
tigate the properties of these processes. To the best of out knowledge all of these methods 
need a scaling relation in order to measure the information included in the series. Another 
well-known method is the level crossing method (LC) in which no scaling feature is explicitly 
required [l-6] and this is the main advantage of this method in estimating the statistical 
information of the series. 

What is the reason that level crossing method has been created? This method was in- 
vented to study the series with different insight. By level crossing, we can measure the 
memory, non-Gaussianity and waiting time (lengt h) ( an average time (length) interval that 
we should wait for an event to take place again |7Hll|). In order to find the way that corre- 
lations enters the level crossing method, we have considered some fractional Gaussian noises 
(fGns) (generalization of ordinary discrete white Gaussian noise) whose Hurst exponents 
and distribution functions are known. Indeed, the so-called Hurst exponent, iJ, gives a 
quantitative measure of the long-term persistence of a signal. In particular, the exponents 
< H < 0.5 and 0.5 < H < 1 correspond to negative (anti-correlation) and positive corre- 
lation, respectively, while H = 0.5 represents an uncorrected Gaussian process. Since fGns 
are well known examples, their comparison with empirical data can be used as a criterion to 
better understand the results obtained from the level crossing method applied to unknown 
empirical data. 

Here, an integrated quantity namely, A^^^^ which represents the total number of upcross- 
ings of a typical series reflects how memory plays role. For a better investigation of the 
memory effects, we have calculated shuffled counterparts of each of underlaying time series 
and compared their associated total number of so-called crossings, (A's/j), with that of given 
by their original time series, N^^., to obtain the percentage of the change in the system. By 
the shuffling procedure, autocorrelations are destroyed. Finally, we have considered the data 
from some markets and mentioned technique has been applied to daily log-returns of three 
time series of financial markets named S'&P500, Dow Jones and Tehran Stock Exchange 
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FIG. 1: Schematic of upcrossing for an arbitrary level, x(t) — x = a. 



(TSE) in the same time interval from 4 Jan 2006 to 4 Jan 2010 [l2l]. Although it is said 
that S&PSOO and Dow Jones have Hurst exponent of about 0.5, the results show that they 
are not exactly white noises. Actually memory exists which have not been observed. Also, 
the source for this inconsistency can be obtained (table [T]). 

This paper is organized as follows. In Section II, we will give a brief history and expla- 
nations concerning the level crossing method. Then, data description and analysis based on 
this method for different fractional Gaussian noises and also its application to some empirical 
data are given in Section III. In Section IV, we present our conclusion. 



II. LEVEL CROSSING ANALYSIS 

For the sake of clarity, we begin with a summary of the LC analysis jl-6|. Consider 
a typical time series {x{t)} and denotes the number of positive difference crossings 
(upcrossings) at level x{t) — x = a in time interval T (see Fig. [1]). The mean value of 
for all time intervals be N^{T) j^, 

KiT) = {ntin, (1) 
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FIG. 2: A schematic of joint probability density function of underlaying series and its derivative 
with respect to dynamical parameter that here is time. The shaded area corresponds to desired 
region satisfying in Eq. (3) 



here (.), represents ensemble average. For a homogeneous (stationary) process, the average 
number of upcrossings is proportional to the time interval T. Subsequently, 



iv+(r) = z/+r. 



(2) 



in which is the average frequency of upcrossings at the level equates io y = x{t) — x = 
a. The frequency parameter i/+ could be deduced from the underlying joint probability 
distributions of y = x{t) — x and y' = = namely p{y,y') Im]. To this 

end, consider a time scale At of a typical sample function, if x{t) — x < a at time t and 
x{t) — X > a at t + At OT alternatively the changes in x{t) is positive in the time interval At, 
there will be a positive crossing of x{t) — x = a. In an other word, two following necessary 
and sufficient conditions should be satisfied to have an upcrossing in time interval. At 



x(t) — X < a 
A [x{t) -x] a 



At 



> 



[xit) 



(3) 
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FIG. 3: Time series of three fractional Gaussian noises with Hurst exponents H, 0.2, 0.5 and 0.8 
from up to down. 

In order to examine whether the above conditions are satisfied at any arbitrary time t, we 
should find how the values oiy = x{t) —x and corresponding time derivative are distributed 
by considering their joint probability density p{y,y')- Suppose that the level y = a and 
interval At are specified. Then, we are only interested in values of y < a and values of 
y' > corresponding the region between the lines y < a and y' > in the plane of 
{y,y') (see Fig. [2]). Therefore, the probability of positive slope crossing of y = a in At is 

oo a 

Y,^y' Yl ^ypiy^y')- (4) 

a-y'At 

We can replace y with a when At — t- 



p{y,y') = p{y = a,y'). (5) 

Since at large values of y and y' the probability density function approaches zero fast enough 
(see Fig. [2]), consequently expression (jl]) may be written as 

a a 
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dy' j dyp{y = a, y') (6) 

'0 Ja-y'dt 

in which the integrand is no longer as a function of y so that the first integral is just 



£.y>dtdyp{lL=a,y' 
is equal to |2|, yj 



p{y = a, y')y'dt. Finally the probability of upcrossing oi y = a 



p{a,y')y'dy' 



(7) 



in which the term p{a,y') is the joint probability density p{y,y') evaluated at y = a. 

The integration of the Eq. ([7]) over all levels demonstrates another quantity, N^^, which 
shows the total number of upcrossings for the data. 



tot 



(8) 



To study the effects of correlations or memory, A^^"^ is evaluated which gives the total number 
of upcrossings of time series when it is shuffled. Here, random permutation is used for 
shuffling the data. The autocorrelations are destroyed by the shuffling procedure. Hence, 
by comparing A^^^^ of the original data with that of computed for shuffled data set, A^^, we 
can obtain the magnitude of correlations in the time series and this gives useful information 
about the time series. In order to study the effect of distribution function in the data with 
non-Gaussian distribution function, the so-called surrogate procedure can be applied. In the 
surrogate method, the discrete Fourier transform (DFT) of the observed time series data 
is computed and then the phases of each complex amplitude of the DFT are replaced with 
independently distributed artificial uniform (— vr, +7r) variates IJ, llSj- The altered DFT is 
then inverse Fourier transformed to generate a surrogate time series. The correlations in 
the surrogate series could be kept unchanged (for more detail see appendix (B) reference 
but the probability function changes to a Gaussian distribution 11 . 
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m\m- To 



obtain information about the effect of the phase randomization procedure on the PDF, one 



15j. For the 



can also check the results of this procedure on the magnitude and sign series 
surrogate time series, iV^ can be calculated which is the total number of upcrossings of this 
time series. 
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FIG. 4: Left panel) Positive level crossings (upcrossings) of different fractional Gaussian noises. 
The average values of all fGns have been shifted to zero and their standard deviations are one. 
Right panel) N^^^ of different fGns with respect to Hurst exponent. Symbols are given by numerical 
analysis and fitting function is given by Eq. ([9]). 

III. DISCUSSION AND RESULTS 

By using the PDF and the correlation function of the data, a stationary series can be 
pictured. Fractional Gaussian noises which are generalization of ordinary discrete white 
Gaussian noise are characterized by their Hurst exponent. Hurst exponent, H, gives a 
quantitative measure of the long-term persistence of a signal. In particular, the exponents 
< H < 0.5 and 0.5 < if < 1 correspond to negative (anti-correlation) and positive 
correlation, respectively, while H = 0.5 corresponds to an uncorrelated Gaussian process. 
For clarifying the level crossing results of unknown empirical data, their comparison with 
the results of well known fGns is emphasized. Fractional Gaussian noises can be generated 
using different methods 20M28[| . Here, we have used the Fourier filtering method which was 



fully described in Refs. [20H22[. In this method, the Fourier components of an uncorrelated 
sequence of random numbers are filtered by a suitable power-law filter in order to introduce 
correlations among the variables. In Fig. El we have plotted schematically three fractional 
Gaussian noises which are from up to down anticorrelated, random and correlated Gaussian 
noises. For a better comparison, their plots have been shifted vertically. 

Level crossing function, is calculated for different fractional Gaussian noises, according 
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FIG. 5: Schematic comparison of positive level crossings for three fGns and their shuffled. 

to Eq. [71 Fig. (jl]) (left panel) shows as a function of level a for different values of Hurst 
exponent. By increasing the Hurst exponent of the fractional Gaussian noises the total 
number of upcrossings decreases (see Fig. (jl]) (right panel)). The behavior of the N^^ 
elucidates the fact that by increasing the Hurst exponent the fluctuation in the time series 
decreases. The solid line shown in Fig. (jl]) (right panel) corresponds to a rational function 
of the following form that has been fitted to the curve of this figure with the goodness of fit 
equals to r = 0.999. 

, 0.668 - 0.615iJ 

JV+ = (Q) 

*°* l- 0.726i/ + 0.0121/2 • V J 

Now it is interesting to separate the correlation and the PDF effects by shuffling and 
surrogate procedures. These procedures are used in the level crossing method to elucidate 
its sensitivity to memory and PDF Since the PDF of the fractional Gaussian noise 

is the normal distribution, there is no more information in the surrogate method. By using 
the shuffling procedure, we can evaluate the memory of the series. In the case of fractional 
Gaussian noises, the total upcrossings (downcrossings) of the shuffled time series is increased 
for the correlated time series {H > 0.5) while this behavior is reversed for the anticorrelated 
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FIG. 6: The relative difference, TZ, as a function of H. Symbols are given by numerical analysis 
and the solid line corresponds to the fitting function which is given by Eq. (jlOp . 

ones (see Fig. ([5])). For the white noise, there is no change in the total positive level crossings 
because there is no correlation in the time series. For a better comparison, their plots have 
been shifted horizontally. By comparing the difference between N^^ and (after shuffling), 
the memory of the time series can be determined. Smaller relative difference denotes that 
the time series is less correlated (anticorrelated). As it is seen (Fig. ([5])) under shuffling 
procedure the total positive level crossings N^^, is increased (decreased) and this confirms 
that the underlying data set is correlated (anticorrelated). In order to quantify the value of 
memory (correlation and anticorrelation) embedded in the date set, we define the relative 
change of the total number of upcrossings by using the A^^^^ for original and shuffled series as 
TZ = \N^^ — N^ot\l^tot- Fig- ©, Tl has been plotted for various values of Hurst exponents 
for different Hurst exponents. Two typical functions have been fitted to the curve of this 
figure which are 




< iJ < 0.5 



0.5 < < 1 



(10) 
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TABLE I: Total positive level crossings of the log-returns of three time series of financial markets 
(5'&P500, Dow Jones and TSE) for the original, shuffled and surrogated time series, and the 
relative difference of the shuffled and the original series one with its own pdf and the other when 
the pdf has been changed to a Gaussian. For better comparison these parameters are also presented 
for a white noise. 



Index 


H 






N+ 

su 


\Nth-Kt\lKt 


\Ku-Kt\lKt 




0.45 


0.49 


0.51 


0.57 


0.03 


0.16 


Dow Jones 


0.53 


0.40 


0.43 


0.57 


0.08 


0.43 


TSE 


0.71 


0.29 


0.49 


0.46 


0.66 


0.59 


White noise 


0.5 


0.56 


0.56 


0.56 


0.00 


0.00 



with the goodness of fit r = 0.999 for each of them. 

The probability density function of stationary empirical data set can be deviated from 
the Gaussian shape, consequently the total upcrossings of the so-called surrogate data is 
a relevant quantity to make our analysis complemented. In order to demonstrate, how TZ 
and the other defined parameters works, we applied this method to log-returns of three time 
series of financial markets which are S'&P500, Dow Jones and TSE. The results are reported 
in table [B For better comparison, we have presented the results for a white noise. Here, we 
have used both the shuffling and surrogate procedures. As seen in the table, for shuffling 
procedure there is not that much change in the total positive level crossing for the two first 
markets which are considered as efflcient markets. But, the shuffling procedure has strong 
effect on TSE which is an inefflcient market. When the surrogate procedure is applied, the 
total positive level crossing for S&PSOO is not that much affected but for Dow Jones and 
TSE this is not the case. This is because their probability density functions are deviated 
from the Gaussian function. In order for them to be comparable to the fGns, they should be 
first surrogated. In fact, correlation and PDF both affect the level crossing method. Here, 
by using both of these procedures we want to consider both their effects. 
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IV. CONCLUSION 



Level crossing method with no required scahng feature is a powerful method in charac- 
terizing the time series. This method not only introduces roughness in a new concept but 
also contains information about waiting time (the time interval to observe an event again, 
statistically). It can also measure the memory and non-Gaussianity. Besides all these in- 
ferences, there was lack of a suitable measure for better understanding the level crossing 
results of unknown empirical data. Fractional Gaussian noises with normal distributions 
which can be indicated by only a Hurst exponent are known samples that could play the 
role of this criterion. They can fill this vacancy by comparing their level crossing results 
with the results of the unknown empirical Data. In this article, the level crossing results of 
fGns are investigated. The deviation of the empirical data results form the fGn with the 
same Hurst exponent is a measure of correlation and non-Gaussianity in the empirical one. 
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